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in the Monte Carlo simulations. 
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1. Introduction 

We have developed a method to derive the (approximate) quark contribution to the fermion 
free energy of QCD on a lattice, at finite temperature and chemical potential, with Kogut-Susskind 
fermions in the flavor basis. This result has been obtained at the lowest order of the nilpotency 
expansion [1-3], an approach which will be outlined in the following. 

For the understanding of our result it is only necessary to anticipate that in our approximation 
the QCD vacuum is dominated by static chromomagnetic fields. At vanishing temperature the 
expression of the partition function we have found is 

2? « J [dU] exp (-S magnetic {U) - S matter {U)) (1.1) 

where S magn etic is the gluon action restricted to spatial plaquettes whose values do not change with 
time, S mat t er is the quark action which depends on the spatial link variables denoted by U and more 
precisely is given by 

= -y£ 2 f i + e(N i -2wah f i)\a^e- 2 > 1 1 + \ (n? + ^4N? + N^J j (1.2) 

where Lq is temporal size of the lattice and is, therefore, the inverse temperature, and should be 
sent to infinity, ji is the chemical potential, d is the step function and N, are the eigenvalues of ./V 
which is (twice) the Dirac Hamiltonian of Kogut-Susskind fermions 



N = -2 Y(i 




j j j j 



(1.3) 



In the definition of this Hamiltonian and are Dirac and taste matrices, and Pj^\v^ are 
projection operators and covariant derivatives 

Pf ] = l -{\®\±Y^®t 5 ^) , V ;. +) = i(^7} (+) -l) , v5.-) = I(l-7; H £/;) (1.4) 

where 7^ ± - > are forward / backward translation operators of one block used to define the Kogut- 
Susskind fermions in the flavor basis and Uj is the j'-th component of U, the spatial link variables 
associated to the blocks of size twice the lattice spacing. 

We notice some features of our expression of the partition function: 

1) In the leading order of our expansion the well known "sign problem" does not arise and 
the quark contribution to the action can be used as a statistical weight in the Monte Carlo 
procedure. 

2) Even though the quark-determinant appears to be non-local, its evaluation requires only the 
knowledge of the eigenvalues of the local Hamiltonian N/2. In the formal continuum limit 
the quark contribution to the action becomes proportional to \N\/2 — }l and we can expect 
that when we approach the continuum non-local effects vanish. 

3) The spatial link variables do not depend on time and the temporal one does not appear at 
all, so that the dimensionality of the system is effectively reduced by one. 
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As in the functional integral formulation the fermionic fields are represented by Grassmann 
variables in a Berezin integral, it is quite difficult, in this formalism, to understand what are the sta- 
tistically relevant configurations which drive the interesting physical phenomena. We have there- 
fore been induced to go back to the operator formulation in the Fock space representation, where 
we used Bogoliubov transformations generating Cooper pairs, whose structure functions can be 
studied in a variational approach, and only afterwards we went back to the functional formulation. 

In general standard Bogoliubov transformations change the terms of an action which are in- 
dividually invariant with respect to some symmetry into terms which are no longer invariant, even 
though the total action remains of course symmetric. It is therefore dangerous to perform ap- 
proximations on the transformed action. To avoid this difficulty in a series of papers [1-3] we 
generalized Bogoliubov transformations in the following way. We perform at each time slice an 
independent transformation whose coefficients are functions of the spatial link variables and of ad- 
ditional bosonic compensating fields. These fields become dynamical fields which describe bosonic 
composites of quarks and antiquarks. The resulting action can be studied in a nilpotency expansion 
whose asymptotic parameter is the index of nilpotency of the composites, which is the number of 
fermionic components in their structure functions. The lowest order is the saddle point approxima- 
tion, in which composites and quasiparticles move in a background field which is the solution of 
the saddle-point equations. 

In the application to QCD at finite chemical potential we assume that the most important 
quark-quark correlations are pairwise, because diquarks are thought to be stable substructures at 
low baryon density (nucleons [4-6], multiquark mesons [7]) and basic constituents of the color 
superconducting phases at high baryon density [8-11]. If this were true diquarks should give con- 
tributions of the same sign to the free energy, and large cancellations among fermion determinants 
at nonzero chemical potential would be due to highly fluctuating, energetically unstable fermionic 
configurations. We then constructed the QCD ground state in terms of diquarks by means of ap- 
propriate (time dependent) Bogoliubov transformations. 



2. First Bogoliubov transformation 

We start from the transfer-matrix formulation of the partition function for Kogut-Susskind 
fermions [12] 

J f, = |[JC/]exp[-5 G (C/)]Tr F | °j\ (f t t ^exp(2 J un B )f m ) | (2.1) 

where 

t t = exp ( vN t u) , V t = exp (u ! In Uqj u + v 1 " In Uq t v) (2.2) 

and Sq is the full gluon action, the baryon number operator, Tr F is the trace on the fermion 
Fock space, u,v are fermion-antifermion canonical annihilation operators, Uqj are temporal link 



variables and N t = N(U t ) is defined in Eq. (1.3) in which the spatial link variables are defined at 
time t. 
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We evaluate the trace on states obtained at each time from fermion coherent states by the 
Bogoliubov transformation 

i 

a = R2(w-^ t v t ) , $ = (v + tf&*)R (2.3) 
R=(l + Jft^-i ^ R=(l+^j^ t )- 1 . (2.4) 

The quasiparticle vacuum has the form of a condensate of fermion pairs \&) = exp (u^^v^) |0). 
Such transformation must be carried out independently at each time slice if we want that the quasi- 
particle operators a, /3 have definite transformation properties. Moreover at each time the matrix 
& must depend on the link variables at this time and in general on compensating fields which 
describe dynamical bosons: 

(^kx 2 =L^(M)(<^x(^)) • (2-5) 

K X 1: X 2 

<p*-(x, t) are bosonic fields with quantum numbers K, and ( &K,x(U t ) ) their structure functions 

\ ' ) Xl,X2 

which must depend on the spatial link variables. Since the Bolgoliubov transformation is unitary 
we could perform a transformation with an arbitrary ^-matrix, that is arbitrary <ps:'s, and we can 
integrate over them in the partition function with an arbitrary measure djj,((p* ,(p). The trace over 
the transformed states in the partition function can be performed exactly yielding its functional 
form 

3T = J [dU]exp[-S G (U)} J JAi(<P>)exp(-%) (2.6) 

where 

s eff = s mesons (j?) - £ «; (v, - j%) a t+1 - A+i (v, - ) A* + /W 2J) « f + afy^p; . 

(2.7) 

We do not have the space to report the explicit expressions of S mesons {^), of the quasiparticle 
hamiltonians ,3% and jf? t and of the coefficients ,J? t '' ' and ' '. 



3. Nilpotency expansion and background field 

Let us consider the composite operator 4>J^ = u*$> K v^ . It is characterized by the index of 

nilpotency, which we denote by 0.% , denned as largest integer such that I <£> K I ^ 0. £1% (usually 
much greater than the number of the internal degrees of freedom of the fermions) is the maxi- 
mum number of composites we can put in the state K. A necessary condition to approximate 
by a canonical bosonic operator is therefore that 0.% >> 1. So the index of nilpotency of com- 
posites which approximate physical bosons can be assumed as asymptotic parameter to set up an 
expansion in its inverse, the nilpotency expansion. To construct such an expansion we look for the 
minimum of S e ff with respect to & neglecting the quasiparticle contribution. This is the saddle 
point approximation which provides the background field J£~. We find the remarkable result that 
jzK 1 - 2 )^) = J^ 2 ' 1 ) ( J£") = 0, meaning that in the background field there is no direct quasiparticle 
-antiquasiparticle mixing. We then set & t = & + 8^ t and expand S mesons in powers of the fluctua- 
tions SJ? t . This results to be an expansion in the inverse of Q. that we call nilpotency expansion [1]. 
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The background field determines the vacuum energy and therefore the phases of the theory while 
the fluctuations 8JP t represent meson fields. 

We found [3] an exact expression of the background field at zero temperature and chemical 
potential, which requires stationarity of gauge fields in the sense that spatial plaquettes are constant 
in time while spatial-temporal plaquettes vanish. We showed that in the saddle point approximation 
color is confined in the quasiparticle sector because quasiparticles propagate only in point-like color 
singlets. We have checked that the nilpotency expansion reproduces correctly the results of a four- 
fermion model both at zero and nonzero chemical potential. 



4. Second Bogoliubov transformation and diquark action 

In the Hamiltonian formalism diquarks are constructed in terms of positive energy states, 
which correspond to quasiparticles in the formalism of the transfer matrix at the saddle point. 
Therefore, at fixed gauge configuration, we construct diquarks as Cooper pairs of quasiparticles by 
a second Bogoliubov transformation 

d = A (&-#tftt) ( 4 _!) 

where 

1 (4.2) 



The matrix @ in an anti-symmetric matrix with the quantum numbers of the diquark field. The 
vacuum of the new quasiparticle operators 6 is 

= exp Q&t^ft^ \W) = exp Q&^ft^ exp (V ^V) |0) , (4.3) 

namely a condensate of Cooper pairs of quasiparticles living in the background 

The saddle point equations for the background field are not changed by the presence of di- 
quarks. Using their solution the diquark field action can be written 

Smatter = y tr {hi (l - le^^) + ^\*{\+&9>) 

1 + 9 {e 2 ^ - iW) # {e 2 ^ - 2W) T | 



iln 

2 



(4.4) 



Lat us represent our operators in the base in which the Hamiltonian of the quasiparticles is diagonal. 
We will restrict ourselves to the so called simple pairing structure functions for the diquark field 

%j = 8 jic E iic ®i , % = % c (4.5) 

in which any quasiparticle state i is associated to one and only one conjugate state i c . The effective 
quark action takes a minimal contribution for states for which <3l\ vanishes or diverges. We will 
denote by i p the states for which | ^ [ = °°. For given chemical potential this action is minimal if 
the i p are all the states for which 

e^-lWi >1. (4.6) 



Introducing the expression of the quasiparticle Hamiltonian [2] we get Eq. (1.2). 

Of course, integrating over the gauge fields with their statistical weight will smooth out the 
distribution of the values of ^, 's. 
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5. Perturbative expansion in the gauge coupling constant 

Suppose that for sufficiently high values of the chemical potential an expansion with respect 
to the gauge coupling constant can be justified 

e^-lW^ l+A + gB + g 2 C. (5.1) 

Assuming simple pairing we get the standard expression of the diquark action 

Smatter ~ - ^ £ { -2fi - In (l +A U + gB u + g 2 Cjj) + Pi (A fi + gB u + g 2 Co) 



2 

i 



-^PiBijPjBji + \ ( ¥ * A t + A* Wi ) | (5.2) 



where 

and ^ 

A,- = - g 2 e iic £ e kkc B ik B ickc yr k (5.4) 
z k 

is the celebrated gap function. By variation with respect to & a gap equation is obtained of the 
standard form, compatible with standard results also in the sense that the gap is dominated by 
chromomagnetic fields with static propagator [13]. 



6. Conclusion 



We have investigated QCD at finite chemical potential guided by the theoretical indications 
that two quarks correlations are important at all baryon densities. We introduced such correlated 
pairs in the formalism of the transfer matrix with lattice regularization by means of two indepen- 
dent Bogoliubov transformations at each time slice. Both transformations at each time depend on 
spatial gauge links and compensating fields at that time. This makes it possible to enforce for 
quasiparticles the same symmetry transformations as for quarks. The first transformation produces 
a background field and quasiparticles, the second yields the diquark field in terms of quasiparticles. 

We have formulated a nilpotency expansion for the effective theory, namely an expansion in 
the inverse of the number of fermionic states in the structure functions of the composites, called the 
index of nilpotency. We have studied the effective action in the saddle point approximation of this 
expansion, which is equivalent to a variational calculation, minimizing the free energy with respect 
to a background and a diquark field. According to the solution for the background field the QCD 
vacuum is a dual superconductor (not color superconductor) from which the chromoelectric field is 
totally expelled (perfect dual Meissner effect) and the fermion Fock space contains only point-like 
color singlets. 

We have derived an equation for the minimum of the quark free energy for any given gauge 
field configuration under the assumption of simple pairing. It cannot be evaluated analytically, but 
there is no sign problem in a numerical simulation. 
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If the effective action in the saddle point approximation is expanded in powers of the gauge 
coupling constant, a gap equation is obtained compatible with standard results. The gap is domi- 
nated by static chromomagnetic fields. 

We think, however, that a perturbative expansion in the gauge coupling constant cannot be 
justified, at least in the simple form in which it is usually done. Schematically, at the baryon 
density at which condensation of molecular diquarks is expected, the gauge coupling constant is 
presumably too large. For baryon densities for which the expansion might be justified, on the 
other hand, we expect a BCS ground state, in which dibaryons have a size much bigger than the 
interquark separations, so that the dependence of the diquark structure functions on the spatial 
gauge fields cannot be ignored. 

We hope that our formulation should give a reasonable approximation to the QCD partition 
function for values of chemical potential of the order of the nucleon mass. Increasing the chemical 
potential we should meet chiral symmetry restoration, which we conjecture should be accompanied 
by the vanishing of the background field. If the phase transition is of first order, to determine its 
location we should compare the free energy evaluated in the present paper with that of the chirally 
symmetric phase in which the background field should vanish. But we notice that the latter cannot 
be simply obtained by setting the background field to zero in our equations. In fact in our saddle 
point approximation we disregarded quasiparticles, on the usual, reasonable assumption that they 
are separated from the vacuum by a large gap. If instead the background field vanishes, the particle- 
antiparticle mixing in the action is again active, and if we construct the diquark field in terms of 
particles, we have no reason to expect that a gap exists for antiparticles. We must therefore proceed 
in a different way that we will illustrate in a future work. 
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